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Abstract 



Based on the system-reservoir description we propose a simple solvable microscopic 

O \ model for a nonequilibrium bath. This captures the essential features of a nonsta- 

O 



tionary quantum Markov process. We establish an appropriate generalization of the 
fluctuation-dissipation relation pertaining to this process and explore the essential 
modifications of the Bloch equations to reveal the nonexponential decay of the Bloch 

►^ ^ vector components and transient spectral broadening in resonance fluorescence. We 

OO 

^S| ' discuss a simple experimental scheme to verify the theoretical results. 
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I. INTRODUCTION 

The dynamics of most of the quantum optical phenomena are based on two fundamental processes; 
the coherent interaction between the system (atom/molecule) and the field mode (classical/quantum) 
and the incoherent dissipation of the system. The latter is traditionally modeled in terms of the 
wellknown system-reservoir theory within the appropriate finite temperature quantum statistical 
scheme [0-|^ • Besides thermal reservoirs, the non-thermal reservoirs QJ^ have also been found to be 
important in connection with the development of correlated emission lasers Q and squeezed light 
fields 1^. The essential underlying assumption about the bath, be it thermal or nonthermal, is 
that it is considered to be in a state of equilibrium throughout the process. Very recently a solvable 
microscopic model for a nonequilibrium bath has been proposed @| to explore classically, the influence 
of an initial nonequilibrium excitation in a complex system on the relaxation of a specific quantity 
of interest. In the present paper we extend this treatment to a quantum optical context. Since 
the initial excitation creates a nonequilibrium energy density fluctuation distribution which imparts 
nonstationarity of the bath, it is expected that optical Bloch equations which take into account of 
both the coherent interaction and the relaxation processes within a simplified description of a two- 
level scheme, are likely to be modified by the nonstationarity of the bath 0,0. Based on a quantum 
version of the model we study this essential modification of the optical Bloch equations and explore 
some of the relevant consequences. 

We thus consider a two-level system in contact with a bath which is not in a thermal Boltzmann 
distribution. This nonequilibrium bath is effectively realized in terms of a semi-infinite dimensional 
broad-band reservoir which is subsequently kept in contact with a standard thermal bath which allows 
the nonthermal bath to relax with a characteristic time scale. The important separation of the time 
scales of fluctuations of the nonequilibrium and the thermal bath is that ||^ the former remains 
effectively stationary on the fast correlation of the thermal noise. The model captures the essential 



features of a nonstationary quantum Markov process. The physical situation that has been addressed 
is the following. At t = the excitation is switched on and the intermediate bath is thrown into a 
nonstationary state. We then follow the coherent dynamics of a classical laser-driven near-resonant 
two-level system interrupted by incoherent emissive processes due to nonequilibrium intermediate 
modes after t > to observe the influence of relaxation of these modes on the transient characteristics 
of the system. We show that the decay of the Bloch vector components is nonexponential in character 
so long as the nonstationarity persists. In addition, the nonstationarity of the bath results in time- 
dependence of the diffusion coefficient which manifests itself in the transient resonance fluorescence 
spectra of the two-level system. The underlying physical mechanism of the transient characteristics 
can be understood with the help of a generalized nonequilibrium fluctuation-dissipation relation 
pertaining to this nonstationary quantum Markov process. In the long time limit one, however, 
recovers the standard Bloch equations and the spectral features. 

The outline of the paper is as follows : In Sec. II we discuss the model for nonequilibrium bath and 
the generalization of the fluctuation-dissipation relation corresponding to the nonstationary process. 
The application of the model to explore the modification of the optical Bloch equations and the 
transient spectral characteristics of resonance fiuorescence have been carried out in Sec. III. In Sec. IV 
we propose a simple experimental scheme to verify the theoretical results on nonexponential decay 
and transient broadening effect. The paper is concluded with a summary of the main results. 



II. RELAXATION OF A TWO-LEVEL ATOM IN PRESENCE OF A NONSTATIONARY 

BATH 

To start with we consider a model two-level atom (the system) coupled to a set of relaxing 
modes considered as a semi-infinite dimensional system which effectively constitutes a nonequilibrium 



reservoir. This, in turn, is in contact with a thermally equilibrated reservoir. Both the reservoirs are 
composed of two sets of harmonic oscillators characterized by the frequency sets {ujj} and {^j} for 
the equilibrium and nonequilibrium bath, respectively. The total Hamiltonian is given by 



Ho = -huoaz + hY^ ^jb]bj + hY^ ^/^cj^C/, 
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+ ^Y1 9^.{(^+c^, + cr-cl) + hY,Yl (^j^ib]cf, + 6jcJ,) . (1) 

The Hamiltonian is essentially a simpler quantum version of the model used in with two-level atom 
as the system. The first term on the right-hand side describes the system mode with characteristic 
frequency Uq. The second and the third terms represent the thermal and the nonequilibrium linear 
modes. The next two terms represent the coupling of the nonequilibrium bath with the system mode 
and the thermal bath where the coupling constants are (^^ and ajf^, respectively. In writing down the 
Hamiltonian we have made use of the rotating wave approximation. 

The Heisenberg equations of motion for the system and the reservoir operators are given by 

cx+{t) = iuJo(T+ - i J2 9f^cl^z , (2) 

cr_(t) = -iuJoa_ +iY,9i^c^(^z , (3) 

a, (t) = -2iY, 9i^c^,a+ + 2^ ^ ^/.cj,a_ , (4) 

hjit) = -iujbj - i Y, (^jt^Cf, , (5) 

c^,{t) = -iVt^Cf, - ig^a_ - i J2 Oij,J^j ■ (6) 

3 

Making use of the formal integral of Eq.(5) for hj{t) in Eq.(6) we obtain 

j 

Jtn 



Taking into consideration that the interference time of Y,j C(jfj_aj^e~'^'^^^*^* ^ is much smaller than the 
time over which the significant phase and amplitude modulation of the linear modes c^{t) takes place, 
the last term in Eq.(7) can be identified as a relaxation term in the usual way m with the damping 
constant 

%, = n a,f,{Q,) a,,{Q,) D{Q,) , (8) 

where D{Q^) represents the density of states of the equilibrium modes evaluated at Q^,. Thus one 
can write down the Langevin equation for the relaxing mode c^ as follows; 

c^(t) = -tn^c^it) - tg^a^t) -Y.llMt) + /mW • (9) 

V 

Here the last term /^(t) represents the usual noise operator arising out of the coupling of the relaxing 
modes with the thermal bath modes as given by 

/,(t) = -^^a,,e--^(*-*«)6,(to) , (10) 

i 

where the reservoir average of f^ify is zero, i.e., 

kfMB = . (11) 

Here by average {0(t))B of an operator 0(t) we mean 

{0{t))B = Tr{0{t)pB} ,wheTe pB = U, exp{{-hu,b]b,)/KT}/[l-exp{-^)] . 

Pb is the thermal operator for initial density matrix for the thermal bath ( initial density for the 
thermal bath {bj}, nonequilibrium bath {c^} and the system are assumed to be factorizable ). 

We now make the following approximations. The cross-terms which involve rapidly evolving 
imaginary exponentials in the summation among the bath modes in Eq.(9) are neglected with respect 
to the diagonals, slowly evolving terms. This secular approximation is the usual one made in the 



context of master equations for baths [0]. It is wellknown [Q that this approximation is vahd in the 
hmit of weak couphng ( \a^^\ <^ \Qiy\ ) and of a 'flat' bath spectrum for which \a^^\ ~ la^^'i/'l- 

Taking into consideration of the above approximations the Langevin equation for the relaxing 
modes, Eq.(9), can be written in the following form, 

c^(t) = -tn^c^it) - tg^a^t) - 7^^c^(t) + /^(t) . (12) 

To explore the influence of an initial excitation of the intermediate reservoir and its relaxation, 
we now consider the evolution of these linear modes c^ in terms of Eq.(12), which allows a formal 
solution of the following form [^ 

cJt) = <(t) + cjto) e(-^^--^M.)(*-*o) _ ^g f dt' e(-^^--^M.)(*-*') a4t') . (13) 

Jto 

The first term on the right hand side in the absence of the coupling of the system mode represents 
the ( long time ) stationary stochastic solution of the form 

c^(t) = c^ e-i[f^.(*-*o)+<^^] ^ (14) 

where the amplitude c^ (operator) and the phases 0^ (c-number) are assumed to be randomly dis- 
tributed 0. The random distribution of phases and amplitudes in the stationary regime makes 
Eq.(13) an instantaneous solution. The second term on the right hand side in Eq.(13) carries the 
information of relaxation of the c^ modes due to their coupling to the thermal bath and is in the 
form of a typical "memory-type term" [see the discussion later]. The latter is not to be confused 
with the usual memory term ( or kernel ) commonly arising out of the frequency dependence of 
friction. The third term on the other hand represents the effect of coupling of the system mode to 
the nonequilibrium reservoir. 

We now substitute this solution (13) in Eq.(2) to obtain the equation of motion for the system 
operator in the usual way as 
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(T- 



_(t) = icuoa+{t) - Ta+{t) + Z^a,{t) 



(15) 



with 



r = TT 5-^(1^0) p(^o) 



(16) 



where the p(fi) represents the density of relaxing intermediate oscillator modes. We assume further 
the weak dependence of 7^ on the modes to perform the integration over Q. 

r can be identified as a dissipation constant of the system mode due to the fluctuation of these 
modes. Also note that 



Zt(t) 






-^ > ,9t^ 



<{t) + cUto) e('^^-^MM)(*-*o) 



is the noise operator for the nonequlibrium bath modes with {Z'^{t))iyji = 0. Here by (0(t))jvK we 



mean 



(0(t))jvR = Tr{Oit)p,} , where p, = I[,exp{i~hn,clc,)/KT}/[1 - exp(--^)] . 

Here pc is the initial thermal density operator for the nonequilibrium bath. 

We proceed similarly to obtain the other equations of motion for system operators, o"_ and aj^ as, 



cr_(t) = —iujQa_{t) — Ta_{t) + Z{t)az , 

&^{t) = -2r(l + a,) + 2Z{t)a+{t) + 2Z^(t)a_(t) 



(17) 
(18) 



Introducing the slowly varying operators as 






-.(Jz 



(19) 



we obtain the following Langevin equations, 

^+(t) = -r^+(t) + 2|t(t)4(t) 

S-{t) = -r^_(t) + 2|(t)^,(t) 

where 



(20) 



eit) = -^J:9. c^,\t)e 



cl{to)e 



(21) 



The nonequilibrium generalization of the fluctuation-dissipation relation is now immediately ap- 
parent. Using Eq.(21) we have 



{^Kmt'))NR = T.9l\{cic;)^ne 



j{n^-ujo)(t-t') 



+ (cUto)c^(to))jvR e*(^M— o)(*-*') e^-^^M^o e-r^.^it+t') 



(22) 



We denote the average photon number of the nonequilibrium bath by 



n{n^,to) = {cl{to)Cf,{to))NR , 



(23) 



where to signifies the dependence of average photon number of the nonequilibrium bath on its initial 
state of preparation. Also the steady state average photon number is given by 



ni^,) = («) 



'NR 



After replacing the summation by integration and 7*^ by an average 7 in Eq.(22) we obtain in the 
usual way 



HHm')) 



NR 



r n(cJo) + e-2T(*-*o) rn{uJo, to)] 5{t - t') . 



(24) 



Eq.(24) and {^{t)) = summarizes the essential properties of the stochastic processes due to 
intermediate oscillator bath modes {c^}. It is important to emphasize that the exponential term in 



Eq.(24) [exp{— 27(t — to)}] does not contain time-difference of the two different instants t and t' over 
which the stochastic process is correlated. Thus this exponential term is not to be confused with 
exp{— 27(t — t')} which normally appears as a typical memory term in correlation function as 

(^(t)^(t')) = Constant x exp[-27(t - t')] 

in a non-Markovian stochastic process. The time-difference of the two instants t and t' in the present 
study appear as an argument of a delta function [6{t—t')] in Eq.(24) due to the use of standard broad- 
band reservoirs, t in the exponential function in Eq.(24) is rather a slow time variable, which makes 
the stochastic process due to {c^}-modes nonstationary. In other words the correlation function in 
Eq.(24) is not invariant under time translation. We are thus concerned here with a nonstationary 
quantum Markov process. This consideration is essential for the application of Onsager's regression 
hypothesis for calculation of spectra with two-time correlation function in the present case as shown 
in the next section. 

Rewriting Tn{uJo,to) in Eq.(24) in terms of a deviation from its steady state value Tn^uo) as 

Tn{uJo,to) = D(to) - rn(c^o) , 

we identify a time-dependent diffusion coefficient D(t) in the last equation (24) as 

Dit) = Tnicuo) + [D{to) - Tn{coo)] e'^^^*-*") . (25) 



We thus obtain 



where we denote 



{i\mt'))NR = Tniuo) fl + re-'^^'-'^n 6{t - t') , (26) 



Tn{uJo) = D{oo) 

_ n{u)o,to) 



Eq.(26) is the desired nonequilibrium quantum generalization of the fluctuation-dissipation re- 
lationship. The classical version of the above equation is given in [^]. This relates instantaneous 
fluctuations of the nonequilibrium bath (which itself is undergoing relaxation at a rate 7 due to its 
coupling with the thermal bath ) to the dissipation of the energy of the system mode through F. 
The nonequilibrium nature of the bath is implicit in the initial preparation which creates an initial 
diffusion coefficient -D(to) and also in the exponentially decaying term. 

To check the consistency of the treatment and to allow ourselves a fair comparison with the 
classical treatment we now make the following comments. 

(i) In the steady state limit one recovers the usual fluctuation-dissipation relation for a thermal bath 
at equilibrium. 

(ii) Eq.(26) can also be expressed in terms of energy density fluctuations of the nonequilibrium 
modes. The energy density which is proportional to the power spectrum centered around c^o is given 
by [h = 1] 

O /■+00 _ _ 

47r J-00 

(28) 



-nn{n) + e-2^(*-*°) 



u{n,to) - -fin(fi) 



It is important to note that t is the slow time variable which is well separated from the time-scale 
of thermal noise. The fluctuations of the noise operator ^(t) is now explicitly determined by the 
nonequilibrium state of the bath modes {c^} through its energy density u{Q,t) at each instant of 
time. In other words the instantaneous nonequilibrium energy density distribution of fluctuating 
modes is related to the friction coefficient of these modes on the system degree of freedom through 
a dynamic equilibrium. The classical version of the above equation can be recovered in the high 
temperature limit (n(n) = l/[exp(f2/i^T) — 1] ~ -^) to obtain 



u{Q,t) = -KT + e-2T(*~*o) \u{Q,to) - -KT 



2 
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2 



(29) 



This classical version was discussed earlier [^] in the context of classical kinetics of complex systems. 
Our quantum generalization is more relevant to quantum optical situations as discussed in the next 
section. 



III. MODIFIED BLOCK EQUATIONS AND TRANSIENT RESONANCE 

FLUORESCENCE 

We have discussed above a simple solvable model for a nonstationary quantum Markov process and 
an appropriate generahzation of the fluctuation-dissipation relation pertaining to this process. Two 
immediate consequences are evident. The first one concerns the modification of decay of the Bloch 
vector components in presence of relaxation of the intermediate bath modes. We show here that 
the decay is nonexponential in nature so long as the nonstationarity persists following the sudden 
excitation. The second one centers around the explicit time-dependence of diffusion coefficient due to 
nonstationarity implied in the fiuctuation-dissipation relations (26). The transient noise spectrum of 
the two-level system is therefore expected to bear this signature of time dependence. With this end 
in view we calculate the physical spectrum of the two-level system in contact with the nonequilibrium 
bath driven by a near-resonant classical monochromatic light field. The Hamiltonian of the coupled 
atom-field system reads as follows; 

H = Ho + h [Va+ e-^'^^(*-*o) + Va^ g^'^^^*"*'')] , (30) 

where Hq is given by Eq. (1) and V represents the amplitude of the classical pump field with 
frequency Uc- Proceeding as before, we obtain the Langevin equations for the slowly varying system 
operators as; 
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-(r - i6)S+ + 2^^{t)S, - 2iVS, 

-2T{S, + I) - tVS+ + iVS^ + ^{t)S+ + eHt)5_ 

-(r + i6)S^ + 2^{t)S, + 2iVS, 



where, S{= ujq — uJc) is the detuning and 






C,{t) is the noise operator as given by 



e(t) = ^E^. kW + c.(to) e(-^^''-^)(*-*«) 






JuJc{t-to) 



The noise is characterized by 



{m)NR = 



1 + re 



-27(t-to) 



5(t-t') . 



(31) 



(32) 



(33) 



(34) 



While considering the above equations we emphasize again the separation of time scales 7 ^ F. 

As a next step, we construct the following Bloch equations for one-time averages from the Langevin 
equations (31) [ this requires the calculation of averages like {C,\t)Sz) which include nonstationary 
contribution involving Eq.(34) as shown in Appendix-A ], 



{Sz(t))NR 



-(r - t5){S4t))^n + 2Tn re-^^\S+{t))MR - 2iV{Sz{t))MR , 
-(r + i6){S-{t))NR + 2T{n + l)re-^^'{S-{t))NR + 2tV{Sz{t))NR , 



-2r 



{l-(2n+l)re-2^*}(5,(t))jvR + 



-iV{S+{t))MR + iV{S^{t))r,R. 



(35) 
(36) 

(37) 



12 



The following comments should be made about the Eqs. (35-37) supplemented by Eq.(34) : 

(i) The exponential term in Eq.(34) results in an effective transient modification of decay rates of all 
the Bloch vector components. In the long time limit one, however, recovers the usual decay rates and 
the standard Bloch equations. It is also interesting to note that the two polarization components 
{S+) and {S-^) decay at different rates in contrast to the usual case of equilibrium bath. 

(ii) The nonstationary contributions in the Bloch equations immediately assert that in absence of 
the driving fields {V = 0) the decay of the polarization components is non- exponential in nature. 
This is reminiscent of what has been observed in the relaxation kinetics of classical complex systems 
where the influence of an initial nonequilibrium excitation of others degrees of freedom of a complex 
system on the relaxation of a specific quantity of interest has been explored. 

(iii) Although the noise correlation in Eq.(34) involves an exponentially decaying term, 6{t — t') 
makes the noise instantaneously correlated. This implies that we consider here a broad-band reservoir 
instead of a colored bath. It is important to note that Lewenstein et.al. in a different context 
have considered earlier the atomic decay in presence of a colored reservoir. They have used the 
modified Bloch equations in non-Markovian form (which involves exponentially decaying terms due 
to the finite response time of the reservoir) and shown how the effects of the colored reservoir can be 
inhibited at large driving fields. Thus the origin of the exponential term in Eq.(34) is different. 

We now turn to the second issue, i.e., the calculation of the transient resonance fluorescent spectra. 
Using matrix notation, the above three equations (35-37) for single time expectation values can be 
put in a compact form as; 

^ = M(t) u(t) + f (38) 

where u(t) and f are the column vectors and are given by 
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{S+(t))NR 



U 






and 



^0^ 



V» / 



(39) 



with 



Mftl 



-(r - i6) + 2rfire~27* -2iV 

-iV -2T + 2Tr{2n + l)e-2^* 

2iV 





iV 
-{T + iS) + 2T{n+l)re 



-27* 



(40) 



/ 



Since the calculation of spectra rests on the evaluation of two-time correlation functions of the 
atomic operators it is essential to examine the validity of quantum regression hypothesis in the 
present context. To this end we note the following points. 

The essential statistical properties of the intermediate bath modes are contained in Eq.(34). This 
equation suggests a differential behavior in time dependence of the two terms. First, the exponential 
time dependence is due to the initial preparation at to and subsequent relaxation at any time [t — to 
) of the intermediate bath modes. So the nonstationary nature is implied in this term. On the other 
hand the 6{t — t') term essentially signifies the correlation of intermediate bath fluctuations, ^(t), at 
times t and t' . The presence of 6{t—t') ensures the broad band nature and hence the Markov property 
of the bath. These statistical considerations, therefore, reveal that the dynamics of the two-level atom 
is acted upon by a nonstationary but Markovian stochastic process due to the intermediate oscillators. 
By using the Langevin description of Heisenberg equation of motion. Lax has proved that Markov 
property implies regression theorem as well the converse. The Markov property is defined by the 
requirement that a Langevin force at time t is uncorrelated to any information at earlier time t'. 
The regression hypothesis is a consequence of this requirement. Lax has specifically shown |p that 
Onsager's original statement for an equilibrium situation is valid even for a nonequilibrium situation 
provided the system is Markovian. The validity of regression hypothesis therefore implies that two- 
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time correlation evolves in the same way as one-time expectation value. 

The equation for evolution of two-time correlation functions is then given by 

d 



dr 



v(t2,r) = M(t2,r) v(t2,r) + F(t2 



(41) 



where 



^ {S+{t^ + r)S^{h))NR^ 



v(t2,r) 



/ 



F(t2 







\ 



V 



-V{S-{t2))NR 





(42) 



{S,{t2 + T)S-{t2))NR. 
\{S-{t2 + T)S^{t2))NR] 

The relevant correlation function fi(t2,r) required for calculation of spectra is the first component 
of the vector v(t2, t) and is given by 



Vi{t2,T) = {S+{t2,T)S4t2)) 



NR 



(43) 



We assume that the atom is initially in its ground state. The Eq.(41) is then solved to calculate the 
correlation function (43). The details are given in the Appendix-B. 

At this point it must be emphasized that since we are dealing with a nonstationary situation 
the standard steady state definition of spectrum is not adequate to describe the transient spectral 
features. We therefore resort to a non-steady state spectrum or the so called 'physical spectrum' 
of the emission from the two-level atoms, where the attention is focused on a dynamic evolution of 
the spectrum following an abrupt excitation of the atom and the intermediate oscillator modes. The 
main reason for studying the time-dependent spectrum is that the familiar power spectrum which 
results from Weiner-Khintchine theorem is not applicable to nonstationary processes. Eberly and 
Wodkiewicz [|l^ have shown that the suitably normalized counting rate of a photo detector can 
be used to define a time-dependent spectrum. This definition allows the influence of the spectrum 
analyzer ( basically a Fabry-Perot interferometer, for example ) to be exhibited in the spectrum so 
that the band limit of the measuring device is appropriately incorporated which makes the spectrum 
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free from ambiguities. It has also been emphasized [10] that when the instrumental width, W, is 
narrow enough such that W <^r the spectrum appears to be qualitatively close to Weiner-Khintchine 



spectrum. This transient spectrum has been used previously on several occasions |]TT|. Following 



Eberly and Wodkiewicz |10| we define the time-dependent spectrum in terms of the time correlation 



function vi(t2,T) as follows; 

S{t, u, W) = 2W Re t dt2 e-'^^*-*^) /*'*' dr e(^-*^Ni(t2, r) . (44) 

JO JQ 



Here t is the elapsed time after the system and the intermediate oscillator modes have been subjected 
to the initial excitation at t = to(= 0), VT is the full width of the transmission peak of the interfer- 
ometer and A(= uj — uo^ is the detuning, or frequency offset of the Fabry-Perot line center above 
the frequency of the field Uc- It is important to note that the time- dependent spectrum is expressed 
in terms of the two integrals. The first integral is over the correlation time r and is actually the 
counterpart of Weiner-Khintchine spectrum band limited by the width W of the measuring device, 
while the second one over ^2 takes into account of the nonstationarity which makes the spectrum 
i-dependent. 

Making use of Eq.(43) in Eq.(44), performing the integration over r and ^2 and extracting the real 
part, we obtain numerically the time-dependent spectrum as discussed below. 

Since the excitation at t = 0, prepares an initial nonequilibrium energy density of the intermediate 
oscillator modes which differs from its equilibrium value, the initial diffusion coefficient -D(O) deviates 
from its stationary long time value D(oo). This deviation is measured in terms of r (= ^,^\ — Ij 
[see Eq.(27)] or equivalently in terms of the ratio of the photon numbers "i^^' ?' . Another quantity of 
interest is the rate of relaxation 7 of the nonequilibrium intermediate oscillator modes due to their 
coupling to the thermal modes. Both r and 7 contribute significantly to the nonequilibrium version 
of fluctuation-dissipation relation [ Eq.(26)] which is essential for understanding the influence of a 
nonequilibrium bath on the transient fluorescence spectrum. In Fig.(l) we plot the physical spectra 
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at three different times after the initial excitation at t = for the parameter set r = 0.4, n = 0.1, 
r = 1.0, 7 = 0.1 under resonance condition 6 = and instrumental linewidth W = 4.0 for a low value 
of field strength V = 2.5. The choice of parameter space is guided by the early work of Eberly and 
Wodkiewicz |T^ on the physical spectra of resonance fluorescence. For a higher value of V{= 10.0) 
the expected three-peak Mollow structure appears as shown in Fig. (2) ( all other parameters kept 
fixed ). The steady state is clearly achieved around t = 10.0. It is interesting to note that in the 
spectra showed in Figs.(l) and (2), the asymptotic condition corresponds to the usual stationary bath 
case, presenting the narrowest line. The transient nature of broadening of the bands is due to the 
creation of a nonequilibrium bath mode energy density through r. Since r is a measure of the excess 
energy gained by the intermediate modes into which the system has to release its energy, the system 
also draws some energy from these modes by virtue of the fluctuation-dissipation relation. The latter 
relation illustrates a dynamical balance of inward flow of energy due to fluctuations from the reservoir 
into the system and the outward flow of energy from the system to the reservoir due to dissipation 
of the system mode. This nonstationary diffusion of fluctuations from the intermediate bath modes 
into the system leads to transient spectral broadening. This persists so long as the nonstationarity 
remains. If the relaxation of these bath modes approaches the time scale of ^ ( with increase of 7 ) 
the broadening effect no longer appears. In the next section we describe an experimental scheme to 
show how this can be realized within the purview of a simple cavity QED experiment. 



IV. DISCUSSION OF AN EXPERIMENTAL SCHEME AND CONCLUSIONS 

Based on a microscopic model for a nonequilibrium bath we have constructed the modified Bloch 
equations which incorporate the effect of nonstationary relaxation and calculated the transient res- 
onance fluorescence spectra of a two-level system driven by a near resonant strong classical fleld. 
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We now discuss a specific system presenting the transient broadening effect studied in tliis paper. 
It is wellknown tliat tlie spectrum of tlie radiation emitted by a strongly driven system is considerably 
modified if the atoms are confined in a cavity. For our purpose the dynamics may be conveniently 
described if one considers a two-level Rydberg atom as a system contained in a cavity ( whose modes 
serve as the intermediate oscillator modes of the present model ). The cavity in turn is weakly 
coupled to the vacuum modes playing the role of the equilibrium reservoir through the cavity losses. 
By sudden sweeping of the resonance of the cavity it is possible to dump an appreciable amount of 
energy on the cavity modes by changing the number of photons abruptly. ( The tuning of cavity in 
studying the emission of strongly driven two-level systems like Ba atoms into the modes of the cavity 



had been experimentally carried out both under adiabatic [|1^ and non-adiabatic |ir3| , p!^ conditions in 
cavity QED experiments [pTl-[13[| . ) This corresponds to the initial preparation of the nonstationary 
state of the cavity modes by changing r = "w^' ?' , see Eqs.(26) and (27) in such a way that the 
energy of these modes becomes suddenly higher than the average energy. Once this nonstationarity 
is attained, the effect of relaxation of the cavity modes on the emission of the strongly driven ( 
externally ) two-level atoms can be monitored by observing the transient fluorescence spectrum. 
Since the atom-cavity interaction ( F, say ~ 100 MHz ) is strong compared to the decay rate of the 
cavity modes ( 7, say ~ 20 MHz ), the separation of time scales as required can be conveniently 
maintained. We also expect to observe the nonexponential decay of emission of the excited two-level 
atom into the modes of an optical cavity so long as the nonstationarity persists. 

We now summarize the main conclusions of this study: 

(i) The microscopic model proposed here may serve as a simple solvable paradigm for a nonstationary 
quantum Markov process. 

(ii) We establish an appropriate generalization of the fluctuation-dissipation relation and its classical 
correspondence pertaining to the above-mentioned process. 

(iii) The origin of nonstationarity (or nonequilibrium nature of the bath) lies in the creation and 
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subsequent relaxation of an energy density fluctuation distribution function of the intermediate bath 
modes following a sudden excitation. 

(iv) Keeping in view of the systematic separation of time-scales involved in the overall dynamics we 
have shown that the decay of the polarization and population inversion components of the Bloch 
vector is non-exponential so long as the nonstationarity persists. 

(v) The nonstationarity of the bath results in time dependence of the diffusion coefficient which show 
up in the transient broadening of the physical spectra of resonance fluorescence. 

(vi) We have outlined a simple experimental scheme within a cavity QED setup to verify the aspects 
nonexponential decay and transient broadening of emissions from a strongly driven two-level system 
in a cavity. 

Since the underlying model of relaxation employed here bears its origin in complex coupled sys- 
tems one may also envisage guest-host systems embedded in a lattice ( where the immediate local 
neighborhood of the guest comprises intermediate oscillator modes and the lattice plays the role 
of a thermal bath ) as typical candidates for experimental realization of such transient fluorescent 
processes. We thus expect the model to be relevant in the context of single molecule spectroscopy 
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APPENDIX A: CALCULATION OF THE AVERAGE {CHT)Sz{T))nr 

To calculate {^'^(t)Sz{t))NR we proceed as follows : 
We have from Eq.(36) 



eit) = -^T.9,\ ct^t) + ct (to)e(''^-^)(*-*o) 



-iujc{t-to) 



This can be written as 

e^(t) = r^(t)+e^^(t) , (Ai) 

where, i^\t) = —iY^^i g/icl} {t)e^^'^''^^^^°^ represents the stationary ( long time ) fluctuation and 

denotes the fluctuations due to the coupling of the system with the relaxing modes. It is essential 
to note that because of the relaxation ^^^(t) noise is nonstationary. It is important to note that the 
separation of time scales of ^'^^(t) and ^^^(t). ^'^^(t) is much faster compared to ^^^(t) and represents 
a Gaussian white noise. Also ^'^^(t) and ^^^(t) are assumed uncorrelated P]. 

Thus we note 

(e^(t) S,{t))^n = {e\t) SMNR+ie'^it) SMnr . 

Since, ^''^(t) is much faster and describes a stationary process, we write 

iCKt) Sz{t))j,n = (r^(t))7v« {SMnr = . 

Thus we have 

i^it) SM)nr = (e^^W S,{t))^n . (A2) 
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Because of the exponential term exp[— 7(t — tg)] in the expression of .^'''^(t), expression (A2) describes 
a nonstationary average which cannot be equated to zero as shown below. 

Making use of the identity 

S,it) = S,it-At)+ [' dt' S,{t') 

Jt-At 

and the expression for Sz(t), we get 

S,{t) = S,{t - At) 

+ f dt' [-2T{S,{t') + \}- iVS+{t') + iVS.{t') + mS+{t') + iKt')S.{t')] 

Jt-At Z 

= S,{t-At) 

+ f dt' [-2T{S^{t') + \}- iVS+{t') + iVS.{t') + e{t')S+{t') + e{t')S^{t') 

Jt-At I 

^i'\t!)S^{tl)^e'\tl)S^{t!)\ . 
We then calculate the average 



{e\t)smNR = f dt' {e\t)eit')s^{t')) 



t-At 



NR 



where we make use of the fact that ^*^(t) and ^"^^{t) are uncorrelated and the operator Sz{t') at 
time t' is not affected by fluctuation at a latter time t. Following Bourret p6| . [T7[] and van Kampen 
]T8| we now make decoupling approximation ( which implies that the correlation of fluctuations C,^(t) 
is much short compared to the coarsed-grained time scale over which the average (5'+) evolves in 
time ) to obtain 

{e\t)e{t')S^it'))^n = {e\t)eit'))^n {S+{t'))NR . 
We then use the fluctuation-dissipation relation for the ^^(t) [ Eq.(37) ] to obtain 

{^Kt) S.it))^n = Tnr e''^' (5+(t));vi? • (A3) 
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We thus take note in passing that identification of a nonstationary part C,'^^{t) (which is not 
invariant under time translation) in Eq.(Al) leads us to a non zero average like (A3). The other 
nonstationary averages are similarly calculated. 

APPENDIX B: CALCULATION OF TWO-TIME CORRELATION FUNCTIONS 



The matrix M as defined in Eq.(40) may be rewritten as 

M(t) = Mo + 2r r e^^T* Mi 



:Bn 



where 



/ 



Mo 



[T-i6) -2iV 
-iV -2T 
2iV 





iV 
{T + 15) 



:b21 



Ml 



n 

2n + 1 
n+1 



The solution of Eq.(41) in terms of Mo and Mi is 



v(t2,r) = exp 



Mqt - -r e~^^(*^+")Mi 



7 



X < v(t2,0) + dt exp 



-Mot + -re~^^(*2+") Ml 
7 



We assume that the atom is initially in its ground state. Then, 

V(t2,0) = 



2V- I 2{Sz)nr) 

— 2{S-)nR 



\ 



^o^ 



, u(0) 



/ 



F(t2 



\' / 



(B3) 



(B4) 



22 



where v(t2,T) and u(t) are defined by Eqs.(42) and (39) respectively. 
Defining a matrix T and a vector g as 

/ 



1 




M 

2 


-| 


, g = 










.°J 



we can write 



v(t2,0) = Tu(t2) 



(B5) 



Hence from Eq.(B4), using Eq.(B5) we get 



v(t2, r) = exp 



MoT - -r e"2^(*2+")Mi 



7 



X [Tu(t2 



dt exp 



-Mot + -re"2^(*2+") Ml 



7 



F(t2 



(B6) 



The solution of Eq.(38) is 



u(t2) = exp 



Mots - -r e-^^^'Mi 



7 



X < u(0) + dt exp 



-Mot + -re-2'^* Ml 

7 



(B7) 



We again define a matrix 



K 



'000^ 

-r 

yO y 



to write 
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F(t2) = K U(t2) . 



(B8) 



Using Eq.(B7) and Eq.(B8) we get from Eq.(B6) , the solution for two-time correlation function in 
terms of the initial condition as follows; 



v(t2, t) = exp Imqt - -re~2T(*2+^)Mi J T J exp 



MnU - -re 



7 



-27*2^,") I u(0) 



+ exp Mot2 re 

V 7 



-27*2 M^j I'' ^^ exp[-|Mc 



r 



t - -re-^^'yii \ f 



7 



+ exp [Mot - -re-2T(*2+")Mi j | f" dtexp [ - JMot - -re^^^^^^+'^Mi l j I 

K iexp (Mot2 - -re-^^^'MA u(0) 

+ exp(Mot2--re"^^*^MiJ /*' rft exp ( - | Mq t - -re'^^^Mi I J f I . (B9) 



24 



REFERENCES 

[1] W. H. Louisell, Quantum Statistical Properties of Radiation (Wiley, New York 1973) ; M. Lax, 
Phys. Rev. 145, 110 (1966) ; M. Lax and H. Yuen, Phys. Rev. 172, 362 (1968) ; G. S. Agarwal, 
Phys. Rev. A2, 2038 (1970) ; R. Graham and H. Haken, Z. Phys. 235, 166 (1970). 

[2] A. Caldeira and A. J. Leggett, Ann. Phys. 149, 374 (1983). 

[3] See for example, a review, G. Gangopadhyay and D. S. Ray, in Advances in Multiphoton Pro- 
cesses and Spectroscopy, edited by S. H. Lin, A. A. Villayes and F. Fujimura (World Scientific, 
Singapore, 1993) vol.8 ; G. Gangopadhyay and D. S. Ray, Phys.Rev. A46, 1507 (1992); A43, 
6424 (1991) ; J. Chem. Phys. 96, 4693 (1992). 

[4] N. Lu and J. Bergou, Phys. Rev. A40, 237 (1989) . 

[5] C. Gardiner and M. J. Collet, Phys. Rev. A31, 3761 (1985); A. E. Ekert and P. L. Knight, Phys. 
Rev. A42, 487 (1990); Z. Ficek and B. C. Sanders, J. Phys. B27, 809 (1994) ; A. S. Parkins, 
Phys. Rev. A42, 4352 (1990). 

[6] J. Ray Chaudhuri, G. Gangopadhyay and D. S. Ray, J. Chem. Phys. 109, 5565 (1998) ; M. 
Millonas and C. Ray, Phys. Rev. Letts. 75 1110 (1995). 

[7] R. Landauer, J. Stat. Phys. 9 351 (1973); 11 525 (1974); 13 1 (1975); D. L. Stein, R. Doering, 
R. G. Palmer, J. L. van Hemmen and R. M. McLaughlin, Phys. Letts. A136, 353 (1989). 

[8] M. Lewenstein, T. W. Mossberg and R. J. Glauber, Phys. Rev. Letts. 59, 775 (1987). 

[9] M. Lax, Phys. Rev. 157, 213 (1967). 
[10] J. H. Eberly and K. Wodkiewicz, J. Opt. Soc. Am. 67, 1252 (1977). 
[11] B. Deb and D. S. Ray, Phys. Rev. A49, 5015 (1994); J. H. Eberly, C. V. Kunasz and K. 



25 



Wodkiewicz, J. Phys. B13, 217 (1980); B. Renand, R. M. Whitley and C. R. Stroud, J. Phys. 
B9, L19 (1977). 

[12] D. Holm and M. Sargent III, Phys. Rev. A33, 4001 (1986). 

[13] Y. Zhu, A. Lezama and T. W. Mossberg, Phys. Rev. A39, 2268 (1989). 

[14] G. S. Agarwal, Phys. Rev. A41, 2886 (1990). 

[15] Ph. Tamarat, B. Lounis, J. Bernard, M. Orrit, S. Kummer, R. Kettner, S. Mais and T. Basche, 
Phys. Rev. Letts. 75, 1514 (1995) ; B. Lounis, F. Jelezko and M. Orrit, Phys. Rev. Letts. 78, 
3673 (1997). 

[16] R. C. Bourret, Can. J. Phys. 40, 782 (1962). 

[17] A. Brissaud and U. Frisch, J. Math. Phys. 15, 524 (1974). 

[18] N. G. van Kampen, Phys. Rep. 24, 171 (1976). 



26 



Figure Captions 



1. Fig.l : Time-dependent resonance fluorescence spectra of the two-level system for different 
dimensionless times with W=4.0, 6 = 0.0, n = 0.1, r=0.4, 7 = 0.1 and V=2.5 (scales arbitrary). 

2. Fig. 2 : Time-dependent resonance fluorescence spectra of the two-level system for different di- 
mensionless times with W=4.0, 6 = 0.0, n = 0.1, r=0.4, 7 = 0.1 and V=10.0 (scales arbitrary). 
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